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RESEARCH  OBJECTIVES 

The  propagation  rate  of  a  premixed  turbulent  flame  is  determined  by  the  complex 
interaction  between  the  incident  turbulent  flow  field  and  the  turbulent  flame  front.  The 
nature  of  this  interaction  and  the  relevant  flow  field  parameters  are  illustrated  schematically 
in  Figure  1.  The  primary  consequence  of  the  turbulence-flame  interactions  is  an  increase  in 
the  flame  area  and  a  proportional  increase  in  the  burning  rate,  due  to  wiinkling  of  the 
laminar  flame  front.  The  local  flame  geometry,  however,  also  affects  the  burning  rate  through 
variations  in  the  local  laminar  flame  speed  caused  by  flame  stretch  effects.  In  addition  to  the 
direct  effect  of  the  incident  turbulence  on  the  burning  rate  through  changes  in  the  flame 
structure,  the  local  velocity  field  produced  by  the  flame  itself  alters  the  turbulence  both 
ahead  of  and  behind  the  flame,  which  in  turn  effects  the  flame  structure  and  therefore  the 
burning  rate.  Thus,  a  complete  characterization  and  understanding  of  turbulence-flame 
interactions  is  essential  for  the  formulation  and  validation  of  phenomenologically  correct 
models  of  premixed  turbulent  flame  propagation. 

The  objective  of  this  research  is  to  develop  an  improved  understanding  of  turbulence- 
flame  interactions  and  their  effect  on  premixed  turbulent  flame  propagation.  More 
specifically,  the  objectives  of  this  research  are  to  experimentally  characterize  and,  in  turn, 
develop  an  improved  understanding  of  the  relationships  between  the  incident  turbulent  flow 
field,  the  upstream  and  downstream  generated  turbulence,  the  turbulent  flame  structure,  and 
the  mass  burning  and  flame  propagation  rates  in  premixed  turbulent  flames.  The  experiments 
are  conducted  in  a  turbulent  flow  reactor  over  a  range  of  pressures,  temperatures  and 
turbulence  conditions  which  include  those  encountered  in  practical  combustion  systems,  such  as 
gas  turbine  combustors.  The  mean  velocity  and  turbulence  are  characterized  both  upstream 
and  downstream  of  the  propagating  flame  using  laser  Doppler  velocimetry.  This  includes 
measurements  of  the  turbulence  intensity,  time  scale,  length  scale  and  energy  spectrum,  as 
well  as  Reynolds  stress  and  vorticity.  The  turbulent  flame  structure  is  characterized  using 
laser  planar  imaging  techniques  and  processed  using  a  fractal  analysis. 

Through  interactions  with  other  researchers,  the  results  of  this  study  will  be  used  to 
develop  improved  and  phenomenologically  correct  models  of  premixed  turbulent  flame 
propagation. 
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FIGURE  1.  Turbulence  -  Flame  Interactions 
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During  the  past  year  experiments  have  been  conducted  in  the  pulsed-flame  flow  reactor 
shown  schematically  in  Figure  2  for  the  case  of  atmospheric  pressure,  unconfined,  propagating, 
one-dimensional  premixed  turbulent  flames.  In  particular,  stoichiometric  propane-air  flames 
have  been  studied  over  a  range  of  turbulence  Reynolds  numbers.  The  measurements  which 
have  been  made  include  ensemble  averaged  velocity,  turbulence  intensity  and  integral  time 
scale  as  a  function  of  time  through  the  propagating  flame.  The  integral  length  scale  has  also 
been  measured  in  the  cold  flow  upstream  of  the  flame.  In  addition,  two-dimensional  flame 
structure  measurements  have  been  made  and  the  resultant  flame  boundaries  processed  using  a 
fractal  analysis.  The  results  of  these  measurements  are  presented  and  discussed  below. 

Measurements  were  made  at  two  turbulence  conditions  which  are  defined  in  Table  1. 

Case  1: 

Case  2: 

The  integral  length  scale,  L,  was  obtained  from  the  two-point  spatial  correlation  measurements 
shown  in  Figure  3,  where  the  integral  scale  is  equal  to  the  area  under  the  spatial  correlation 
curve.  A  value  of  46  cm/sec  was  used  for  the  laminar  flame  speed,  s^,  in  stoichiometric 
propane-air. 

A  set  of  twelve  two-dimensional  flame  structure  measurements  were  obtained  for  each 
turbulence  condition.  The  individual  images  were  processed  to  define  the  flame  boundary 
separating  the  burned  and  unburned  gases.  The  flame  boundary  was  then  processed  using  a 
fractal  analysis  procedure  suggested  by  Mandelbrot  [1]  and  Sreenivasan  [2]  in  order  to  obtain 
the  fractal  dimension.  A  typical  two-dimensional  flame  boundary  and  the  corresponding  fractal 
analysis  are  shown  in  Figure  4.  For  all  of  the  measurements,  the  fractal  analysis  yielded  a 
linear  curve.  This  indicates  that  the  flame  structures  are  self-similar  which  is  a  characteristic 


Table  1 


n’/st  =  0.26 

L  =  3.7  mm 

u’/s/  =  0.54 

L  =  4.3  mm 

of  fractal  geometries.  For  this  particular  method  of  fractal  analysis  the  resultant  fractal 
curve  is  defined  by  the  following  relationship, 

A  *  2  «3‘D  L 

where  A  is  the  area  defined  by  all  points  which  lie  within  a  distance  e  from  the  flame 
boundary,  L  is  the  length  of  the  flame  boundary  for  a  given  f,  and  D  is  the  fractal  dimension. 
Therefore,  the  slope  of  a  plot  of  log  A  versus  log  e  is  equal  to  3-D,  which  can  then  be  used 
to  determine  the  fractal  dimension. 

The  fractal  dimension  was  determined  for  each  measurement  as  described  above  and  then 
an  average  fractal  dimension  was  calculated  for  each  turbulence  condition.  This  average 
fractal  dimension  is  plotted  versus  an  integral  length  scale  based  turbulence  Reynolds  number 
in  Figure  5.  Also  shown  are  the  fractal  dimensions  which  were  obtained  previously  for  other 
turbulence  conditions.  The  combined  results  show  a  monotonic  increase  in  fractal  dimension 
with  increasing  turbulence  Reynolds  number,  where  the  line  at  D  =  1.37  represents  a  high 
Reynolds  number  limit  for  the  fractal  dimension  as  proposed  by  a  number  of  researchers  [3,4], 
These  results  are  significant  in  that  they  clearly  demonstrate,  even  at  low  turbulence  Reynolds 
numbers,  the  fractal  nature  of  premixed  turbulent  flame  structures.  They  also  suggest  the 
existence  of  a  functional  dependence  of  fractal  dimension  on  turbulence  Reynolds  number. 

Such  a  relationship  would  be  important  for  developing  a  model  of  turbulent  flame  propagation 
which  explicitly  accounts  for  changes  in  flame  structure. 

The  LDV  velocity  measurements  are  made  at  a  fixed  location  as  a  function  of  time  as  the 
turbulent  flame  passes  through  the  measurement  location.  Measurements  are  typically  made 
over  400  or  more  flame  events  and  the  individual  velocity  time  records  are  shifted  to  match 
flame  arrival  times.  The  shifted  velocity-time  records  are  then  used  to  calculate  the  ensemble 
averaged  velocity,  turbulence  intensity  and  temporal  autocorrelation  function  as  a  function  of 
time  through  the  propagating  flame. 

The  ensemble  averaged  velocity  results  are  efcown  in  Figure  6.  Both  cases  shew  the 
expected  sudden  increase  in  velocity  after  flame  arrival,  however,  the  velocity  increase  is  less 


than  ''ne  would  expect  due  to  the  unconfined  nature  of  the  flame.  Also,  note  that  both 


flames  show  a  slight  decrease  in  the  velocity  ahead  of  the  flame,  where  again  this  can  be 
attributed  to  the  fact  that  this  is  an  unconfined  flame.  That  the  velocity  decrease  appears  to 
start  much  earlier  or  further  ahead  of  the  higher  turbulence  intensity  case  flame  is  primarily 
due  to  the  difference  in  the  laboratory  flame  of  reference  propagation  rates  of  the  two 
flames.  The  higher  turbulence  intensity  case  flame  propagated  much  more  slowly  which  serves 
to  "stretch"  to  time  axis  on  these  plots.  Once  the  propagation  velocity  is  accurately 
measured,  it  will  be  possible  to  convert  this  data  from  temporal  to  spatial  profiles. 

The  ensemble  averaged  turbulence  intensity  results  are  shown  in  Figure  7.  Again,  the 
higher  turbulence  intensity  case  result  can  be  thought  of  as  having  an  expanded  time  axis  due 
to  its  slower  propagation  velocity.  Both  flames  show  clear  evidence  of  an  increase  in 
turbulence  intensity  ahead  of  the  flame.  This  is  more  clearly  resolved  in  the  higher 
turbulence  intensity  case,  where  nearly  a  factor  of  two  increase  is  observed.  Flame  arrival  is 
indicated  by  the  sudden  increase  in  turbulence  intensity.  With  respect  to  the  upstream 
turbulence  intensity,  there  is  a  tenfold  and  a  fivefold  increase  in  the  low  and  high  turbulence 
intensity  cases,  respectively.  It  is  perhaps  more  relevant  to  note  that  the  absolute  change  in 
turbulence  intensity  is  approximately  the  same,  i.e.,  1  m/sec,  in  both  cases.  This  is  consistent 
with  the  limited  amount  of  similar  data  which  is  available  from  other  experiments  [5,6]. 

The  velocity-time  data  is  also  used  to  calculate  the  ensemble  averaged  temporal 
autocorrelation  coefficient.  A  typical  temporal  autocorrelation  curve  is  shown  in  Figure  8, 
where  the  area  under  the  autocorrelation  curve  is  equal  to  the  integral  time  scale.  In  Figure 
9,  the  integral  time  scale  is  plotted  versus  time  through  the  propagating  flame  for  both 
turbulence  cases.  In  both  cases,  there  is  a  decrease  in  the  integral  time  scale  immediately 
ahead  of  the  flame,  while  behind  the  flame  the  time  scale  is  observed  to  increase  in  one  case 
and  decrease  in  the  other.  Whether  or  not  these  results,  particularly  behind  the  flame,  are 
meaningful  will  require  additional  measurements.  Assuming  that  Taylor’s  hypothesis  is  valid  in 
this  flow,  one  can  calculate  the  integral  length  scale  from  the  integral  time  scale.  The 
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results  obtained  ahead  of  the  flame  are  shown  in  Figure  10.  If  these  are  compared  to  the 
length  scale  measurements  obtained  from  the  two-part  spatial  correlation  measurements,  one 
finds  that  they  differ  by  a  factor  of  approximately  three,  which  is  not  unexpected  since  the 
isotropic,  homogeneous,  equilibrium  turbulence  assumptions  required  for  Taylor’s  hypothesis  are 
not  valid  in  this  flow. 
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